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Abstract 

This talk gives a very short overview of some of the important physical 
phenomena observed at the Relativistic Heavy Ion Collider (RHIC). 
The emphasis is put on the multiplicities, hard probes and the proper- 
ties of the initial state. 

1 Introduction 

The Relativistic Heavy Ion Collider (RHIC) is a unique machine designed to create a very high 
energy density over an extended region as a result of nuclei collisions. This process enables to 
investigate the collective phenomena in Quantum Chromodynamics. In particular, one hopes to 
create the quark-gluon plasma which is the state of deconfined quarks and gluons. According to 
cosmology such state existed at the very early Universe. RHIC machine is capable of accelerat- 
ing and colliding different hadronic systems: proton-proton, nucleus-nucleus, deutron-nucleus at 
a wide range of centre-of-mass energies, from 19.6 to 200 GeV per nucleon pair. Just after the 
collision, a very high quantity of quarks and gluons is created. It is estimated that, the formation 
time for the initial density is about 0.35 fm/c. After that, the system thermalizes very quickly 
and reaches thermodynamic equilibrium. Estimates suggest that this happens rather quickly, 
with very short thermalization times of the order of 1 fm/c. The assumption of thermalization is 
vital for the application of the hydrodynamics [1] which is used to describe the expansion of the 
system. When the system expands and cools down, the quarks and gluons form hadrons which 
finally reach the detectors. The phase transition from the quark-gluon plasma to hadrons oc- 
curs at rather large value of strong coupling which is beyond the applicability of the perturbative 
methods. Nevertheless, one can explore thermodynamic properties of QCD using lattice meth- 
ods, see for example [2]. In Fig. [T] we show the result of the lattice calculations [2] for the energy 
density divided by T 4 as a function of temperature. The results clearly show the phase transition 
at the critical temperature T c of about 173 MeV. The critical energy density corresponding to 
this temperature is about e c ~ 0.7 GeV/fm 3 . The energy density at RHIC [3] can be calculated 
from the transverse energy density at midrapidity via Bjorken formula: (e) ~ ^iff > where 
A is the overlap area for the colliding nuclei. The average energy density depends crucially 
on the estimates of the time r at which it is evaluated. For thermalization times in the range 
0.6 — 1.0 fm/c the average energy density is about 9.0 — 5.4 GeV/fm 3 . This is well above the 
critical density obtained from lattice calculations, compare Fig. [I] It is interesting to note that, the 
energy density of the cold nuclear matter is about 0.15 GeV/fm 3 . In Fig. [I] the energy densities 
probed by RHIC and LHC are also indicated. Notably, the calculations signal that already at SPS 
energies, the transition from hadron phase to quark-gluon plasma occured. Lattice calculations 
enable to probe the phase diagram of QCD, shown schematically in Fig. [2] The vertical axis is 
the temperature T and the horizontal one is the baryo-chemical potential. 
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Fig. 1: Energy density /T 4 as a function of the temperature. The blue arrow with e/T 4 is the Stefan-Boltzman limit. 
Figure by F. Karsch [4]. 



The region above the dashed-solid lines is the quark-gluon plasma phase, whereas the 
region below, at small temperatures and baryo-chemical potential, is the hadron phase. The 
transition between the two is a smooth crossover. Both RHIC and LHC probe the quark-gluon 
phase and the transition region at small values of the baryon chemical potential. Further to the 
right, for higher values of the baryo-chemical potential the critical endpoint is expected and the 
transition becomes of the 1st order. At high values of the baryo-chemical potential and smaller 
temperatures, new phase appears, the color superconductor. This phase is not accessible at the 
high energy colliders. 
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Fig. 2: Phase diagram of QCD. 
The fact that RHIC reached quark-gluon plasma domain in the nucleus-nucleus collisions 



suggest that, at least in some of the measured observables, one can expect to see dramatic effects 
as compared with proton-proton collisions. In reality the situation is quite complicated. Indeed, 
large differences are seen, for example the suppression of px distributions at large values of the 
transverse momenta with respect to the (scaled) proton-proton collisions. On the other hand, bulk 
properties, like total multiplicities and rapidity distributions, have quite similar shapes and energy 
dependences as compared with the scaled proton-proton measurements. We are going to review 
some of the observations performed at RHIC and discuss the phenomenological descriptions. 



2 Multiplicities 

In Fig. [5] from [3] the measurements of the total multiplicity in AA, pp(pp), e + e~ collisions, 
scaled by the number of participating nucleon pairs iV part /2 are shown as a function of the 
increasing energy. The smooth rise with energy is well described by the In 2 s behavior over wide 
range of the energies. What is striking, is the fact that the data for scaled multiplicity for nucleus- 
nucleus collisions lie on top of the data points for e + e~ collisions. Proton-proton data lie lower, 
most probably due to the leading particle effect. In the ppipp) collisions, lots of the energy is 
taken by the quark spectators into the forward region, and only a fraction of the energy is used for 
the production of the secondary particles. The proton-proton data can be superimposed onto the 
nucleus-nucleus and e + e~ data when the energy is rescaled by a factor of 1/2. This universality 
of the multiplicites indicates that the bulk of the produced particles depends only on the total 
energy (and iVp art ) and not the species of the colliding particles. 
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Fig. 3: Total multiplicity scaled by the number of participating pairs as a function of the energy. Compilation of data 
from e + e~ , pp(pp), AA and dA collisions. Figure by PHOBOS collaboration from [3]. 



The measured multiplicity at mid-(pseudo)rapidity in nucleus-nucleus collisions increases 
slower with energy , as In s. This simple linear in In s extrapolation to the LHC energies predicts 
that N 2 - ^\ v=0 ~ 6 — 7. Various theoretical predictions for LHC energy are larger than this 
simple estimate and span a wide range, up to nearly 40 for the value at midrapidity. This is prob- 
ably connected to the fact that in most of these calculations some semi-perturbative component is 
present which results in a power-like increase of the multiplicity rather than the logarithm. This 
large uncertainty of the phenomenological extrapolations reflects our rather limited theoretical 
understanding of the energy dependence of multiplicities. 



2.1 Extended longitudinal scaling 

PHOBOS collaboration performed measurements of the rapidity distributions for various ener- 
gies (and systems) and found that the distributions exhibit limiting fragmentation property which 
is also called extended longitudinal scaling [5]. This means that when viewed in the rest frame of 
one of the projectiles the (pseudo)rapidity distribution becomes a function of only rf = rj—Y\, earn 
where Ibeam is the rapidity of the beam. Therefore the rapidity distribution in the regime around 
7/ ~ is dominated by the fragments of the broken target whereas the fragments of the projectile 
move with increasing velocity as the energy is further increased (to study these fragments one has 
to go to the rest frame of the projectile rj + Ibeam)- The limiting fragmentation also requires that 
the interaction between the target and the projectile does not depend appreciably on the energy. 
This scaling means that the rapidity distributions must be determined very early in the collision, 
most probably by the initial state. 




Fig. 4: Rapidity distributions for the nucleus-nucleus collisions as a function of the shifted variable r( = r\ — ibeam 
for various center-of-mass energies. Curves are from calculation using the CGC nonlinear equations [6]. 

Fig. |4] illustrates this phenomenon. The curves shown in figure are obtained from the 
calculation [6] based on the kx -factorization approach together with the evolution via nonlinear 
equation for the gluon density. In this approach the rapidity distribution is evaluated as the 
convolution of two parton densities corresponding to the target and the projectile respectively. In 
the fragmentation region, the target parton density is evaluated at rather large values of Bjorken 
x whereas the projectile density at rather small values of x. The essential part of this calculation, 
which enables to reproduce the observed scaling, is the fact that the rapidity distribution in the 
target fragmentation region is dominated by the initial state parton density of the target, probed 
at large values of the Bjorken variable x. At these values the parton density possesses Bjorken 
scaling,i.e. does not depend on the scale at which it is evaluated. The projectile density, which 
depends on rather small values of x and small scales, is saturated since it is evaluated from 



the evolution equation which takes into account nonlinearities important at high partem densities. 
Therefore this density does not depend much on the x values and hence the center of mass energy. 
We see that even though the approach is semi-perturbative it does capture the essential physics 
features necessary to reproduce the limiting fragmentation phenomenon. 



3 Hard probes 

To explore the properties of the quark-gluon plasma created in the high energy nucleus-nucleus 
collisions experiments at RHIC measured the suppression of the production of the high px parti- 
cles and heavy quarks. These are excellent probes of the produced medium, and due to the large 
difference in scales (high px as compared to the bulk of low px particles ), it should be in princi- 
ple possible to employ the perturbative methods. One usually quantifies the effect of the medium 

by evaluating the ratio Raa = (jv ? ii ) d d^ A jdr]dpt • ^ s seen fr° m Fig- [5 from [7] the Raa ratio 
is significantly below 1 even at the very high values of px- This phenomenon is called the jet 
quenching which is interpreted as the interaction of the produced jet (or rather leading high px 
particle) with the produced medium. The interaction results in the significant energy loss of the 
leading high px particle. In fact the effect is so large that it indicates that the high px particles 
which reach the detectors are emitted from a relatively thin outer shell of the high density region. 




Fig. 5: Raa for the produced hadrons as a function of their transverse momentum. Figure from [7]. 



To describe the jet quenching one usually starts from the standard collinear factorization 
formula for the exclusive hadron production in the vacuum 

°£t* h = fA{xi,Q 2 ) ® f B (x 2 , Q 2 ) ® a{xi,xz, Q 2 ) D^ h + 0(A/Q) , (1) 

where Ja, Jb are parton distribution functions evaluated at scale Q 2 , a is the partonic cross 
section and D^~* h is the fragmentation function from parton / to hadron h. To evaluate the 
production in the medium one convolutes the above expression with the so called quenching 
weight [8] 

<r£t h = <r£-+ h ®P(AE,L,q), (2) 



which depends on the length of the medium and the parameter q which acts as a transport coef- 
ficient for the medium. One has to emphasize that although in the proton-proton collisions the 
factorization formula ([I]) is well established, (up to higher twist corrections) it is not proven to 
hold in form ([2]) for collisions which involve large nuclei. Final state interactions could in prin- 
ciple affect the partonic cross section a in a non-factorizable way. It is therefore an assumption, 
that after the partons have been produced in the hard scattering process, one can further factorize 
their interactions into one universal function P. The curves in Fig. [5] are from [8] which is cal- 
culation based on the theoretical framework described above. The quenching parameter q was 
fitted to get the best description of the data. 

The production of the heavy quarks is another excellent process. In principle the calcula- 
tion framework is the same as above. One does expect some differences due to the fact that in the 
vacuum the radiation patter for heavy quarks is different from the light quarks. Due to the their 
mass heavy quarks should radiate less than the light quarks at angles smaller than 6q = m/E 
where m is the mass of the heavy quark and E is its energy. This is called the dead cone ef- 
fect [9]. Therefore the net effect would be less supression (larger Raa) than for the light quarks. 
The measured suppression for the heavy quarks is of the same order as for the light quarks, and 
the theoretical predictions (for example [10]) do not quite predict such a large value. One has to 
emphasize that the heavy quark production in pp (from electrons emerging from semi-leptonic 
decays) is underestimated by the NLO perturbative QCD calculation by a factor of about 2 for 
PHENIX data and by about factor 5 for the STAR data. Therefore the process with the heavy 
quarks calls for better understanding, possible taking into account various effects : collisional 
energy loss, bottom/charm ratio or factorization breaking. 

We conclude our discussion of the hard probes with the description of the new calculational 
methods which employ the AdS/CFT correspondence to evaluate the jet quenching parameter. In 
this approach [11], the expectation value of the Wilson loop, as an average over the medium 
is evaluated using the string/gauge duality in the limit of the strong coupling constant. This 
expectation value is obtained as an exponential of the string action evaluated at the minimum 

(W F (C))=ex V (iS(C)-iS ) , 

where S is the Nambu-Goto action with the metric on the 4 + 1 dimensional AdS space. In the 
case of the finite temperature, the corresponding metric is that of the AdS Schwarzschild black 
hole. In the strong coupling and the multicolor limit the problem becomes classical, i.e. reduces 
to finding the extremum of the action. It is possible, using this method, to evaluate the value of 
the jet quenching parameter which in this case is qsYM = 5 GeV 2 /fm. We note however, that 
all these results can be derived only for the case of the N=4 SYM theory. 

4 Initial state 

The processes described above clearly indicate the presence of the dense medium in the final 
state. The natural question arises whether one can also observe the effects coming from the initial 
state, namely the wave function of the colliding nuclei. BRAHMS collaboration performed a 
measurement of the high pt suppression as a function of rapidity for dA collisions [12]. Whereas 
at midrapidity no suppression is observed for this process, the Rcp clearly shows a decreasing 



trend when moving into forward rapidity. This phenomenon was quite successfully described by 
the models which involve saturation effects in the gluon density or in general by the Color Glass 
Condensate model [13]. In CGC the basic object is the wave function of the nucleus. Obviously 
the complete knowledge of such wave function requires methods which go beyond that known 
in preturbation theory. One can nevertheless calculate the variation of this wave function with 
energy. This is governed by the renormalization group equation which can be derived from the 
Feynman graphs in the leading logarithmic approximation. At very small values of Bjorken 
x one expects the fast growth of the gluon density within the nucleus. CGC model together 
with the renormalization group equations predicts that this growth should be tamed whenever 
the x becomes sufficiently small. The transition between the fast growth and the regime where 
the recombination effects for the gluons become important is governed by the saturation scale Q s 
which is a function of the Bjorken x. Thus the saturation scale provides with a dynamical cutoff at 
low values of x and at low scales. It is the prediction of the CGC model that the Rqp ratio should 
decrease at forward rapidities [14]. CGC model has been successfully used to describe various 
observables in heavy ion collisions: multiplicites, rapidity distributions (mentioned already in the 
previous section) [15] and also Rcp ratio. One has to emphasize though that there are several 
questions concerning strict applicability of this approach to the RHIC data. The values of x are 
not very small for the RHIC kinematics, the formalism correctly incorporates only gluons and 
it has been so far only used at leading order whereas higher orders are known to be very large. 
Nevertheless, the CGC approach, mostly due to its interesting properties in the infrared regime, 
remains a very attractive approach both theoretically and phenomenologically and its predictions 
should be further confronted with the experimental data. 

Conclusions 

We have discussed a selection of the phenomena measured at the Relativistic Heavy Ion collider. 
Clearly, due to our space and time limitations the list presented here is by no means exhaustive. 
The bulk properties as shown in the measurements of the multiplicities are very similar to that 
measured in the simpler systems: proton-proton or even in e + e~ collisions. The extended longi- 
tudinal scaling of rapidity distributions indicate the importance of the initial state. Hard probes 
in form of the high px particles or heavy quarks signify the presence of the strongly interacting 
medium. This is further corroborated by the observation of the strong elliptic flow in peripheral 
collisions. The theoretical descriptions based on hydrodynamics have been quite successful in 
describing the hadron spectra and the anisotropy. Also calculations which employ the perturba- 
tive methods supplemented by the rescattering or recombination(saturation) effects are able to 
describe the bulk of the data at large values of the transverse momenta. Nevertheless, despite 
these incontrovertible successes in phenomenology, the RHIC data still constitute a significant 
challenge for a theory and call for a more complete and coherent description within QCD. 
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